syms n;
y = (1 + n)^(1/n);
limit(y,n,0)
t=[0.01:0.01:100]
z = subs(y,n,t)
plot(t,z)
